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A class of simplified measures is constructed to capture 
the key features of generic spatio-temporally chaotic systems. 
A combined analytical and numerical investigation allows us 
to establish the scaling behaviour of the fractal dimension in 
open systems. Our results improve a previous conjecture and, 
what is more important, furnish a clear framework for both 
numerical and analytical checks of the underlying assump- 
tions. 
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The structure of the invariant measure in generic 
spatio-temporal chaotic systems is still far from being un- 
derstood. The only basic conclusion that is undoubtedly 
accepted within the physicist community is the exten- 
sivity of the fractal dimension, i.e. that the number of 
active degrees of freedom is proportional to the system 
size . The number of degrees of freedom per unit vol- 
ume is the so-called dimension density, a quantity that 
can be derived from the Lyapunov spectrum, through the 
well known Kaplan- Yorke formula. 

However, as soon as finite subsets of (in principle) infi- 
nite systems are considered (i.e., an open-system point of 
view is adopted), it is immediately far from obvious how 
to characterize the probability distribution of generic ob- 
servables. The common belief is that by looking at the 
system with a sufficiently coarse-grained resolution, one 
is not able to distinguish between a closed and an open 
system The stochastic-like action of the external 

world (i.e. the rest of the chain) which activates oth- 
erwise stable degrees can be resolved only if we look at 
the dynamical system with a sufficiently high resolution. 
This was first suggested by Pomeau ||, who postulated 
that the effect of distant (in real space) degrees of free- 
dom on the dynamics in a specific site decreases exponen- 
tially fast with the distance. However, even though this 
statement looks rather plausible, no convincing argument 
has yet been presented, which justifies this hypothesis in 
terms of the actual behaviour of perturbations. 

The same problem has been considered in the spirit 
of time-series analysis In particular, Tsimring || 

formulated a conjecture about the amount of information 
contained in finite samples of temporal signals generated 
by spatio-temporal chaotic regimes. In fact, this problem 



logically descends from the understanding of the struc- 
ture of the invariant measure in open systems, since the 
behavior of a local observable over a finite time can be 
obtained by applying the evolution operator to station- 
ary spatial configurations in a finite region (the so-called 
light cone). For this reason, we shall restrict our discus- 
sion to this latter context in which Korzinov and Rabi- 
novich Bhave proposed a conjecture similar to those in 
in Fiefs . [lJJ^] . Specifically, they claim that the effective 
(i.e. finite scale) dimension D e in subsystems of length L 
(for the sake of simplicity, we consider one-dimensional 
systems) depends on L and on the observational scale e 
as 



vd 2 

D e = dL lne - A, 

V 



(1) 



where d is the dimension density, 77 is the Kolmogorov- 
Sinai entropy density , v is the propagation velocity of 
disturbances, while A is a non-bettcr-spccificd parame- 
ter. Unfortunately, the derivation of the above formula 
is based on several assumptions that cannot be directly 
checked. Moreover, it is rather unlikely that more accu- 
rate numerical simulations will ever provide data clean 
enough to draw definite conclusions. It is therefore com- 
pelling to make some progress on the theoretical side even 
at the expense of introducing strong simplifications. This 
is the route already undertaken in Ref. Q], where the 
limit case of weakly coupled maps has been considered. 
Here, rather than attempting to prove that the invariant 
measure of some dynamical system has a given struc- 
ture, we have preferred to construct a class of measures 
fulfilling the basic requirements for a space-time chaotic 
regime and yet are simple enough to be handled and 
controlled. The key approximation consists in assuming 
that the support of the probability distribution is a linear 
subspace, so that we can use a global approach such as 
singular-value decomposition to get at once information 
on the structure of the invariant measure on all possible 
length scales. Singular-value decomposition as a tool for 
analysing space-time chaos has been already profitably 
applied to experimental || as well as to numerical || 
data, allowing to identify the relevant modes. However, 
the validity of these results is limited by the unavoid- 
able presence of nonlinearities which definitely induce a 
bending of the support as well as local nonuniformities 
in phase-space. Here, by referring to suitable linear sub- 
spaces, we automatically get rid of these effects and can 
use a global methodology to extract local information. 
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We shall consider a scalar process x(i) defined on a 
lattice of length N, with the label i denoting the spatial 
position. A generic configuration is constructed as the 
linear combination of some modes 



x(i) 



D 

E 

k=l 



p k e k (i), 



(2) 



where the p k s are i.i.d. random variables and ek{i) de- 
notes the ith component of the kth mode. In other words, 
the coefficients p k represent the coordinates of the pro- 
cess x(i) in the basis of modes e k . The number of modes 
D = dN can be read as the "fractal" dimension of the 
measure, while d is the dimension density which is as- 
sumed to be independent of the length, i.e. an intrinsic 
characteristic of the underlying "dynamical system" . 

For the sake of simplicity, we first restrict our analysis 
to Fourier modes, 



efc(i) 



1 v Y if fc = l 

^2/Ncos(a k i) if k > 1 and odd (3) 

y/2/Nsm(a k i) if k even, 



where a k = 2n[k/2]/N and [•] denotes the integer part. 

Even if very simple, such a hypothesis is not simplis- 
tic. In fact, let us recall that Fourier modes are the 
stable/unstable directions in a chain of Bernoulli maps 
x t +i(i) = Mod[a(iJWt(i-l) + (l-2(i)xt(i)+tiXt(i+l)), 1]. 
If one chooses the local slope a and the coupling strength 
[a such that the Kaplan- Yorke dimension density JlJ is 
equal to d, the active degrees of freedom correspond pre- 
cisely to the above mentioned D Fourier modes. 

The problem we want to address concerns the struc- 
ture of the projection S N (L,d) of the global invariant 
measure onto the lower-dimensional space corresponding 
to a sub-chain of length L <C N. We expect that the 
resulting distribution extends along all L directions, i.e. 
that its dimension coincides with the space dimension L. 
Nevertheless, we also expect S N (L,d) to be very thin 
along those directions corresponding (in the closed sys- 
tem) to the most contracting directions (i.e., those modes 
e k with k > D Q). 

As we are dealing with a linear subspace, the extension 
of S N (L, d) along the various directions can be obtained 
by means of the standard orthogonal decomposition of 
the correlation matrix 



Cg(L,d) 



(x(i),x(j)) =J2(Pk) e k(i)e k (j), 



(4) 



fc=i 



where (•) denotes an ensemble average and the last equal- 
ity can be easily proved by substituting Eq. (||). It is 
convenient to assume that the p k s are uniformly dis- 
tributed within the interval [— 3 1 / 3 , 3 1 / 3 ] so that (p 2 k ) = 
31/3 

J x 2 dx = 1. The eigenvalues X 2 (1 < I < L) of the 
matrix C (L,d) represent the mean square extensions 
of S N (L,d) along the axes defined by the corresponding 
eigenvectors. 



It is possible to obtain an analytic expression for the 
correlation matrix. By substituting Eq. (Q) into Eq. (Q), 
and assuming that D is odd, we find, after some simple 
manipulations, that 



(D-l)/2 



(•5) 



fc=i 



It is immediately seen that Cu = D/N = d, while 
a general expression for the other entries can be ob- 
tained by using the relation Y^k=i c °s(kx) — cos^(n + 

l)x/2j sia(nx/2)/ sin(x/2) and other simple trigonomet- 
ric formulas, 



sm 



Ci UL,d) 



iVsin 



((i-j)ir/N 



(6) 



In the limit N — > 00 (i.e. for an infinitely extended 
system) but fixed dimension density d and observational 
length L, the correlation matrix converges to 



C«(L,d) 



sinni- j)d7rj 

7T(« - j) 



(7) 



Similar calculations for even D show that Cy converges 
to the same limit, which thus holds in full generality. 

Since the matrix entries depend only on \i— j\, the val- 
ues along all diagonals are constant (Toeplitz structure) 
and the correlation matrix is completely determined by 
its first row or column. Although this structure is rather 
simple, all attempts to diagonalise analytically the cor- 
relation matrix failed. Therefore, we have been obliged 
to perform numerical investigations with high precision 
(IO- 32 ). 

The most convenient way to look at the results is by or- 
dering the eigenvalues Xf from the largest to the smallest 
one. The qualitative behaviour for d = 0.5 and L = 20 
can be seen in Fig. 1, where we have reported the linear 
extensions Xi of the projected measure along the orthogo- 
nal axes represented by the eigenvectors of the correlation 
matrix. We see that approximately 10 values are close to 
1, while the others are almost negligible. This confirms 
that, on a coarse-grained scale, an open system looks like 
a closed system with the same dimension density d. Nev- 
ertheless, the eigenvalues beyond the 10th are not equal 
to zero: this is the result of the coupling with the pseudo- 
random process generated by the external chain. 

From the viewpoint of "fractal" -dimension analysis, we 
can interpret these results by saying that the effective di- 
mension corresponding to a fixed observational resolution 
e = lambdcLi, is D e = i. The expression for the effective 
dimension can thus be obtained by inverting the expres- 
sion for the eigenvalue distribution, 



D e (e) = i(X). 



(8) 
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In the specific case reported in Fig. 1 , the dimension seen 
for e ~ 0.1 is still smaller than 13. Only for e- values as 
tiny as 10~ 7 , the full space dimension (20) is recovered. 

This problem has a meaningful interpretation also in 
the context of linear time-series analysis, since we are 
addressing the question of how the average Fourier spec- 
trum of a stochastic signal (observed in a window of 
length L) converges to the limit shape for increasing L. 
The difference is that here the spectrum of the signals 
are obtained by means of singular value decomposition 
rather than finite Fourier transform. 

Our main goal is to extract the relevant asymptotic 
behaviour of D e on the window-length. As it can be seen 
in Fig. 2, where In Xi/L is plotted versus i/L for different 
choices of L, there is a clear evidence of an asymptotic 
scaling regime such as 

lnXi = -LF(i/L), (9) 

where the function F(x) is identically zero for x < d, 
while it increases monotonously for x > d. From a funda- 
mental point of view, it is important to have a clear idea 
about the "critical" behaviour for x ~ d. In particular, it 
is necessary to understand what a kind of singularity is 
present in that region. The numerical analysis indicates 
that the right derivative F'(x) computed in x = d con- 
verges to a finite value. Accordingly, for x > d, we can 
expand F{x) in a power series, 

oo 

F(x) = J2Pj(x-d) j (10) 

3=1 

By substituting Eq. ( |l0| ) in Eq. (^|) and retaining only 
the first two terms we find the approximate expression 

hxe^-0 x {D e -dL)-^{D e -dL) 2 , (11) 

where we have identified A with e and i with D e . Before 
commenting on the physical implications of the above 
relation, it is important to discuss finite-size corrections. 
A technical analysis exploiting the symmetry properties 
of A 2 suggests that the deviation from the asymptotic 
value of F(x) around i = dL are equal to — (\n2)/2L 

By including such corrections in Eq. ( |11[ ) and thereby 
inverting, we find that 

„ , In e 62 „ „ t In 2 , . 

D e = dL ^-(ne) 2 . (12) 

The structure of this expression reduces to that of the 
conjecture in Ref. |J (i.e. Eq. (|l|) provided that the 
squared logarithmic term is negligible, i.e. provided that 
L 3> lne. One can easily check that this inequality is 
not always satisfied in direct computations of the dimen- 
sion, so that it would be worth re-analyzing numerical 
data in order to obtain reliable estimates of the various 
coefficients. 



A further difference with Eq. (|]) resides in the coeffi- 
cient in front of the logarithmic term: we find 1/ 0\ that 
is to be confronted with vd 2 /r). We are able to determine 
0i only numerically (even though rather accurately) and 
it turns out to be approximately equal to 0.6 indepen- 
dently of d in the interval [0.2,0.8]. The natural model 
to compare our expectations with Eq. ([!]) is again the 
chain of Bernoulli maps in which case, one can deter- 
mine analytically the Lyapunov spectrum pl| and easily 
compute d, rj, and v. Tests made for different choices 
of the slope a and of the diffusive coupling fi show that 
vd 2 jr\ is always significantly larger than l//?i (by almost 
a factor 2). This inequality hints at a possible consis- 
tency of the two approaches. In fact, Eq. ( p^ ) can at least 
be viewed as a lower bound for the effective dimension, 
since it has been obtained by neglecting nonlinearities 
(which enter as discontinuities in Bernoulli maps) which 
can only contribute to increase the dimension on finite 
scales. Anyhow, it would be very important to establish 
whether the discrepancy is to be attributed to a failure 
of the conjecture in Ref. or to ingredient neglected 
in our treatment. In favour of the former hypothesis, we 
claim that, as long as nonlinearities imply only large scale 
phenomena such as bendings and nonuniformities, they 
should not affect the asymptotic behaviour described by 
the scaling function F(x). Indeed, from the point of view 
of an information-theoretic approach (which is implicitely 
that one adopted by deriving an expression for the effec- 
tive dimension), there is no difference between a straight 
and a bended manifold, provided that the bending is not 
too strong |L3|]: the same number of boxes is needed to 
cover the set. 

Aside from the role of nonlinearities, the validity of our 
results could be challenged on the basis of the special 
assumptions made to construct the invariant measure. 
Therefore, we have decided to progressively remove some 
of the limitations. We have started by assuming that the 
average amplitude of the first (nonzero) D Fourier modes 
is not constant but goes continuously to 0. In this case, 
one can again find an explicit expression for the corre- 
lation matrix; numerical studies performed for different 
choices of the Fourier amplitudes indicate a convergence 
towards the same function F{x) but stronger finite-size 
corrections (i.e. larger A). 

A further criticism could be that the Fourier basis is 
a rather special choice, contrasting with the observation 
that Lyapunov vectors are localized E^] . We have there- 
fore considered also bases made of exponentially local- 
ized vectors. This has been done by implicitely referring 
to Bernoulli maps with quenched disorder (i.e. by ran- 
domly fixing the slope of the local maps). In this case, 
the correlation matrix can only be constructed numeri- 
cally and for lattices of finite length N; moreover, there 
is the additional difficulty of the dependence on the re- 
alization of the disorder. Surprisingly enough, we find 
again a reasonable agreement with the results obtained 
for the Fourier basis with the same system-size (provided 
that the geometric average of the eigenvalue spectra is 
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FIG. 1. Extensions along the orthogonal directions for an 
object with dimension density d — 0.5 and an embedding 
dimension L — 20. 

FIG. 2. Scaling of the eigenvalue spectrum of C(L,d) for 
different spatial lengths L but fixed dimension density d — 0.5. 
The increasingly large plateaus seen for large values of L are 
a numerical artifact signaling that the limit of computer ac- 
curacy has been reached. 



taken). Accordingly, we are led to conjecture that the 
behaviour displayed by the Fourier basis is quite univer- 
sal and seemingly independent even of the localization 
properties of the basis. 

Summarizing, we have found a rather general scaling 
law expressing the dependence of the effective dimension 
on the size of an open system and on the observational 
resolution. As the result has been derived under the as- 
sumption of a linear structure for the invariant measure, 
it is now crucial to test its validity in generic models, 
where nonlinearities play certainly an important role. 
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